PLACING FEATHERS ON A SURFACE 

BACKGROUND OF THE INVENTION 
The present invention relates to computer 
generated graphics, in particular, the modeling and 
5 rendering of feathers using a computer. 

Computer graphics are used in many- 
different applications including computer games, 
movies and web pages. Graphically modeling and 
rendering of feathers for a realistic simulation of 
10 birds and other feather objects is a difficult task. 
In general, graphical feathers are generated on an ad 
hoc basis that consumes a significant amount of an 
artist's time to study and draw feathers. Thus, a 
systematic method for rendering feathers is useful. 
15 SUMMARY OF THE INVENTION 

The present invention provides a computer 
implemented method for placing feathers on a surface. 
The method includes providing a surface having a 
plurality of vertices and establishing a growing 
20 direction for each of the plurality of vertices on 
the surface. Feathers are placed on the surface based 
on the plurality of vertices and the growing 
direction. In further embodiments, the feathers can 
be placed using a recursion or interpolation 
25 algorithm. 

Another aspect of the present invention is 
a method for placing feathers on a surface including 
providing a surface having a plurality of vertices. 
Each vertex has a growing direction. The method also 
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includes performing a recursive algorithm to place a 
feather at each vertex. The algorithm includes 
finding a growing direction for vertices in the 
growing direction of the vertex. If the feather at 
5 the vertex collides with another feather, then the 
growing direction of the vertex is adjusted until 
there is no collision. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates a general computing 
1 0 envi ronmen t . 

FIG. 2 illustrates an exemplary feather. 

FIG. 3 illustrates a block diagram of an 
exemplary framework for a system. 

FIG. 4 illustrates an exemplary user 
15 interface for modeling and rendering of individual 
feathers . 

FIG. 5 illustrates example feathers 
rendered using the system illustrated in FIG. 3. 

FIG. 6 is a schematic model of a feather 
20 for sampling a feather texture. 

FIG. 7 illustrates rendering a feather 
according to an embodiment of the present invention. 

FIG. 8 illustrates an exemplary framework 
for placing feathers on a surface. 
25 FIG. 9 illustrates placing feathers on a 

wing and a tail skeleton. 

FIGS. 10A-10D illustrate placing feathers 
on a polygon of a surface . 

FIG. 11 illustrates a first bird and a 
3 0 second bird. 
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DETAILED DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS 

Prior to discussing the present invention 
in greater detail, an embodiment of an illustrative 
environment in which the present invention can be 
5 used will be discussed. FIG. 1 illustrates an example 
of a suitable computing system environment 100 on 
which the invention may be implemented. The computing 
system environment 100 is only one example of a 
suitable computing environment and is not intended to 

10 suggest any limitation as to the scope of use or 
functionality of the invention. Neither should the 
computing environment 100 be interpreted as having 
any dependency or requirement relating to any one or 
combination of components illustrated in the 

15 exemplary operating environment 100. 

The invention is operational with numerous 
other general purpose or special purpose computing 
system environments or configurations. Examples of 
well known computing systems, environments, and/or 

20 configurations that may be suitable for use with the 
invention include, but are not limited to, personal 
computers, server computers, hand-held or laptop 
devices , multiprocessor systems , microprocessor-based 
systems, set top boxes, programmable consumer 

25 electronics, network PCs, minicomputers, mainframe 
computers, distributed computing environments that 
include any of the above systems or devices, and the 
like. 

The invention may be described in the 
30 general context of computer-executable instructions, 



such as program modules, being executed by a 
computer. Generally, program modules include 
routines, programs, objects, components, data 
structures, etc. that perform particular tasks or 
5 implement particular abstract data types. The 
invention may also be practiced in distributed 
computing environments where tasks are performed by 
remote processing devices that are linked through a 
communications network. In a distributed computing 

10 environment, program modules may be located in both 
local and remote computer storage media including 
memory storage devices. 

With reference to FIG. 1, an exemplary 
system for implementing the invention includes a 

15 general purpose computing device in the form of a 
computer 110. Components of computer 110 may 
include, but are not limited to, a processing unit 
120, a system memory 130, and a system bus 121 that 
couples various system components including the 

20 system memory to the processing unit 120. The system 
bus 121 may be any of several types of bus structures 
including a memory bus or memory controller, a 
peripheral bus, and a local bus using any of a 
variety of bus architectures. By way of example, and 

25 not limitation, such architectures include Industry 
Standard Architecture (ISA) bus, Micro Channel 
Architecture (MCA) bus, Enhanced ISA (EISA) bus, 
Video Electronics Standards Association (VESA) local 
bus, and Peripheral Component Interconnect (PCI) bus 

30 also known as Mezzanine bus. 
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Computer 110 typically includes a variety 
of computer readable media. Computer readable media 
can be any available media that can be accessed by 
computer 110 and includes both volatile and 
5 nonvolatile media, removable and non-removable media. 
By way of example, and not limitation, computer 
readable media may comprise computer storage media 
and communication media. Computer storage media 
includes both volatile and nonvolatile, removable and 

10 non- removable media implemented in any method or 
technology for storage of information such as 
computer readable instructions, data structures, 
program modules or other data. Computer storage 
media includes, but is not limited to, RAM, ROM, 

15 EE PROM, flash memory or other memory technology, CD- 
ROM, digital versatile disks (DVD) or other optical 
disk storage, magnetic cassettes, magnetic tape, 
magnetic disk storage or other magnetic storage 
devices, or any other medium which can be used to 

20 store the desired information and which can be 
accessed by computer 110. Communication media 

typically embodies computer readable instructions, 
data structures, program modules or other data in a 
modulated data signal such as a carrier wave or other 

25 transport mechanism and includes any information 
delivery media. The term "modulated data signal" 
means a signal that has one or more of its 
characteristics set or changed in such a manner as to 
encode information in the signal. By way of example, 

30 and not limitation, communication media includes 
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wired media such as a wired network or direct -wired 
connection, and wireless media such as acoustic, RF, 
infrared and other wireless media. Combinations of 
any of the above should also be included within the 
scope of computer readable media. 

The system memory 130 includes computer 
storage media in the form of volatile and/or 
nonvolatile memory such as read only memory (ROM) 131 
and random access memory (RAM) 132. A basic 
input/output system 133 (BIOS), containing the basic 
routines that help to transfer information between 
elements within computer 110, such as during start- 
up, is typically stored in ROM 131. RAM 132 
typically contains data and/or program modules that 
are immediately accessible to and/or presently being 
operated on by processing unit 120. By way of 
example, and not limitation, FIG. 1 illustrates 
operating system 134, application programs 135, other 
program modules 136, and program data 137. 

The computer 110 may also include other 
removable/non-removable volatile/nonvolatile computer 
storage media. By way of example only, FIG. 1 
illustrates a hard disk drive 141 that reads from or 
writes to non-removable, nonvolatile magnetic media, 
a magnetic disk drive 151 that reads from or writes 
to a removable, nonvolatile magnetic disk 152, and an 
optical disk drive 155 that reads from or writes to a 
removable, nonvolatile optical disk 156 such as a CD 
ROM or other optical media. Other removable /non- 
removable, volatile/nonvolatile computer storage 



media that can be used in the exemplary operating 
environment include, but are not limited to, magnetic 
tape cassettes, flash memory cards, digital versatile 
disks, digital video tape, solid state RAM, solid 
state ROM, and the like. The hard disk drive 141 is 
typically connected to the system bus 121 through a 
non-removable memory interface such as interface 140, 
and magnetic disk drive 151 and optical disk drive 
155 are typically connected to the system bus 121 by 
a removable memory interface, such as interface 150. 

The drives and their associated computer 
storage media discussed above and illustrated in FIG. 
1, provide storage of computer readable instructions, 
data structures, program modules and other data for 
the computer 110. In FIG. 1, for example, hard disk 
drive 141 is illustrated as storing operating system 
144, application programs 145, other program modules 
146, and program data 147. Note that these components 
can either be the same as or different from operating 
system 134, application programs 135, other program 
modules 136, and program data 137. Operating system 
144, application programs 145, other program modules 
146, and program data 147 are given different numbers 
here to illustrate that, at a minimum, they are 

different copies. 

A user may enter commands and information 
into the computer 110 through input devices such as a 
keyboard 162, a microphone 163, and a pointing device 
161, such as a mouse, trackball or touch pad. Other 
input devices (not shown) may include a joystick, 
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game pad, satellite dish, scanner, or the like. 
These and other input devices are often connected to 
the processing unit 120 through a user input 
interface 160 that is coupled to the system bus, but 
5 may be connected by other interface and bus 
structures, such as a parallel port, game port or a 
universal serial bus (USB) . A monitor 191 or other 
type of display device is also connected to the 
system bus 121 via an interface, such as a video 

10 interface 190. In addition to the monitor, computers 
may also include other peripheral output devices such 
as speakers 197 and printer 196, which may be 
connected through an output peripheral interface 195. 

The computer 110 may operate in a networked 

15 environment using logical connections to one or more 
remote computers, such as a remote computer 180. The 
remote computer 180 may be a personal computer, a 
hand-held device, a server, a router, a network PC, a 
peer device or other common network node, and 

2 0 typically includes many or all of the elements 
described above relative to the computer 110. The 
logical connections depicted in FIG. 1 include a 
local area network (LAN) 171 and a wide area network 
(WAN) 173, but may also include other networks. Such 

25 networking environments are commonplace in offices, 
enterprise -wide computer networks, intranets and the 
Internet . 

When used in a LAN networking environment, 
the computer 110 is connected to the LAN 171 through 
30 a network interface or adapter 170. When used in a 
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WAN networking environment, the computer 110 
typically includes a modem 172 or other means for 
establishing communications over the WAN 173, such as 
the Internet, The modem 172, which may be internal 
5 or external, may be connected to the system bus 121 
via the user-input interface 160, or other 
appropriate mechanism. In a networked environment, 
program modules depicted relative to the computer 
110, or portions thereof, may be stored in the remote 

10 memory storage device. By way of example, and not 
limitation, FIG. 1 illustrates remote application 
programs 185 as residing on remote computer 180. It 
will be appreciated that the network connections 
shown are exemplary and other means of establishing a 

15 communications link between the computers may be 
used. 

In order to gain a better understanding of 
the present invention, an exemplary feather is 
illustrated in FIG. 2. Feather 200 is defined by a 

20 left outline curve 202 and a right outline curve 204. 
Additionally, feather 200 includes a rachis 206 and a 
plurality of barbs 208, for example a left barb 210 
and a right barb 212. The geometry of feather 200 is 
defined by its rachis and the barbs distributed along 

25 both sides of the rachis. Vanes on both sides of 
rachis 206 carry the plurality of barbs. Each barb 
includes a series of barbules, which are projections 
off of the barbs. 

FIG. 3 illustrates a block diagram of an 

3 0 exemplary framework for a system 250 in accordance 
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with an embodiment of the present invention. In 
system 250, several inputs are provided by a user in 
order for system 250 to model feathers. In order to 
model individual feathers, geometry input 252 and 
5 texture image input 254 are provided to system 250. 
Input 252 includes various parameters defining a 
geometric shape of the feather. The parameters of 
input can 252 include outline curves 256, rachis 
curve 258, barb curves 260 and random seed 262. 

10 Texture image input 254 is used to define the texture 
and color of the feather. Inputs 252 and 254 are 
provided to an interface 264. Interface 264 provides 
an interactive interface for a user to selectively 
alter inputs 252 and 254. Additionally, interface 264 

15 provides inputs 252 and 254 to a feather rendering 
system 266. To render individual feathers, system 266 
includes a bi-directional texture function (BTF) 
module 268 and a parametric L-system module 270. The 
L-system module 270 is coupled to BTF module 268 in 

20 order to output and render synthetic feathers. 

FIG. 4 illustrates an exemplary user 
interface 264 for modeling and rendering of 
individual feathers. The interface includes a rachis 
window 3 02 that defines the rachis curve of the 

25 feather. A barb window 3 04 illustrates barb curves 
with views from a perspective parallel with the 
rachis curve and perpendicular with the rachis curve. 
Window 306 is a sample feather texture image that has 
been input by a user. For example, a user may scan in 

3 0 or otherwise provide a sample feather that is to be 
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modeled and rendered by system 266. Additionally, 
window 308 illustrates outline curves of the feather. 
Each of the rachis curve in window 302 , the barb 
curves in window 304 and outline curves in window 308 
5 may be interactively changed by a user in order to 
control the overall shape of the feather. Once these 
curves have been defined, system 266 is able to 
render individual feathers. 

System 266 uses L-system module 270 to 

10 produce a base structure for the feathers. A feather 
can be regarded as a branching structure composed of 
repeated units called modules. An L-system represents 
the development of a branching structure by 
productions. A production replaces a predecessor 

15 module by several successor modules. A production can 
either be context-free and depend only on the module 
replaced, or be context-sensitive, in which case the 
production depends on the replaced modules as well as 
its immediate neighbor modules. Context-free 

20 productions can be of the form: 

id: pred: cond succ 
where id is the production label, and 
whereas pred, cond and succ are the predecessor, 
condition, and successor, respectively. The 

25 production is carried out only if the condition is 
met . 

Given the rachis, barb, and outline curves 
as entered using interface 264, feather geometry can 
be modeled using a parametric L-system as follows: 



30 
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co : R(O) 

p, : R(i) : i <N - [B L (i,0)] [B R (i,0)]R(i+l) 
p 2 :B L (ij):j<M L ^ B L (i,j + l) 
p 3 :B R (i,j):j<M R - B R (i,j + 1) 
5 (1) 



where N defines the length of the feather 
as well as the density of the barbs at each side of 
the rachis while M L and M R define the lengths of the 

10 left and right barbs respectively. The axiom R(0) 
generates a feather based on the rachis and barb 
curves. Production pi produces a small segment of the 
rachis according to the rachis curve and grows a barb 
on each side of the rachis using recursion. 

15 Production p 2 creates a small segment of the left barb 
according to the left barb curve while production p 3 
proceeds similarly on the right barb with a right 
barb curve. Feather 340 in FIG. 5 shows a feather 
created using the above equation. 

20 However, equation (1) above ignores the 

interaction between neighboring barbs and thus 
feathers appear unnatural. A feather geometry 
generated by equation (1) looks plausible but too 
regular. For a real feather, two lateral sets of 

25 barbs within the vane of the feather interlock the 
feather together. The interlocking is important for 
flight, keeping the air from rushing right through 
the feather. When the interlocking system of a 
feather is disturbed, as when a twig brushes through 



a feather, random gaps form between the barbs on the 
same side of the rachis. 

Neighboring barbs cling to each other by 
little hooks called cilium on the ends of the 
barbules if the total external force exceeds the 
holding capabilities of the cilium. To simulate this 
effect, external forces are introduced into our 
parametric L-system based on random seed parameter 
262 as follows: 

co : R(0, 0, 0) 

pi :R(i,F L ,F R ):i<N&&F L < F 0 && F R < F 0 -> [B L (i,0)] 

[B R (i,0)]R(i+l,F L +F e , F R + F e ) 
p 2 : R(i, F L , F R ) : i < N && F L > F 0 && F R < F 0 - 

[O L B L (i,0)] [B R (i,0)] R(i+1, 0, F R + F e ) 
p 3 :R(i,F L ,F R ):i<N&&F L < F 0 && F R > F 0 - [B L (i,0)] 

[O R B R (i,0)]R(i+l,F L +F e , 0) 
p 4 : R(i, F L , Fr) : i < N && F L > F 0 && F R > F 0 - 

[O L B L (i,0)] [O R B R (i,0)] R(i+1, 0, 0) 
p 5 :B L (ij):j<M L - B L (i, j + 1) 
P6:B R (i,j):j<M R - B R (i,j + l) 

(2) 

where F L and F R are the total external 
forces on the left and right barbs respectively. 0 L 
and 0 R are directional rotations of the left and right 
barbs in response to F L and Fr. The productions pi and 
p 4 produce a portion of the feather that, for each 
step along the rachis curve, F L and F R exceeds a 
threshold force F 0 exerted by the cilium, the left 
(right) barb is rotated by a random angle 0 in a 
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direction determined by F L (F R ) . The rotation of the 
barb is assumed to be within the tangent plane 
defined by the tangent vectors of the rachis and barb 
at the point where the rachis and barb intersect. The 
5 random rotation angle is computed as 0=A,rk(s) where 
rk(s) is the k-th random number generated with random 
seed and can be a user-defined constant. After the 
rotation F L (F R ) starts to accumulate again from zero. 
The random seed of each feather can be saved so that 

10 its shape remains the same every time it is rendered. 
FIG. 5 contains a number of feathers 341-344 created 
with random gaps between the barbs. The different 
feathers 341-344 are rendered by changing various 
feathers geometry parameters. 

15 Equation (2) provides a feather structure 

including random gaps. Adding texture and color to 
the structure provides realistic rendering of the 
feathers. In order to provide texture to the feature 
structure, BTF module 268 captures the mesostructure 

20 and the directional radiance distribution at each 
point on the feather surface. A BTF is a 6D function 
T(x # y, 6 V , <|> v/ 0i, <|> i) , where (0 V/ § v ) is the viewing 
direction v and (0i, (f>i) is the lighting direction 1 
at surface point (x,y) . The BTF defines variation in 

25 pixel intensity based on a number of parameters. 

To calculate the BTF, a geometry model is 
built for the barbs and barbules as shown in FIG. 6 
according to the structure of the barbs and barbules. 
FIG. 6 illustrates a schematic model 360 of three 
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barbs 362 and a plurality of barbules 364 projecting 
from the barbs 362. This structure is rendered for 
all viewing and lighting settings. In one embodiment, 
the rendering is done offline so that a complicated 
5 geometry and sophisticated lighting models can be 
used. This model is opaque with both diffuse and 
specular reflections. 

A sample line 366 of the BTF can be taken 
along an axis perpendicular to the barbs 3 62 to 

10 obtain a 5D BTF T b b(x, 0 V / <K, ©i/ 4>i) for some 
constant y 0 . A 5D BTF suffices for rendering the 
actual 6D BTF of a feather because of the spatial 
arrangement of barbs and barbules. The mesostructure 
of barbs and barbules is rendered such that fine- 

15 scale shadows, occlusions, and specularities are 
well -captured in the BTF T b b- The rendering can be 
done offline by using a ray-tracer. 

The above model of feather mesostructure 
has a number of advantages. First, the off-line BTF 

20 calculations allow capturing of a complicated 
mesostructure and directional radiance distribution 
at each surface point. Second, the BTF can model 
additional effects such as oil-film interference and 
iridescence, which is important for a class of 

25 familiar birds such as hummingbirds and ducks. 
Finally, a level-of -detail rendering can easily be 
supported with a BTF. As a result, the feather 
mesostructure can easily be simulated. 

When rendering a feather, parametric L- 

30 system module 270 generates the feather at run-time. 
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By generating the feather at run-time, the storage 
requirement is modest for each feather because only 
its L-system parameters and the random seeds are 
stored; details such as barb curves (polylines) and 
5 the random gaps on the vane (the feather blade) are 
generated at run-time. After the feather structure is 
rendered using L-system module 270, BTF module 268 is 
used to efficiently draw the barb-barbule 
mesostructure on the feather structure and thus 

10 achieve realistic rendering for a wide range of 
viewing distances. 

FIG. 7 illustrates the rendering of a 
feather. The feather L-system describes a barb 400 by 
a polyline 402 with vertices {x lf x 2 , ... , x n } . After 

15 polyline 402 is generated, texture is added by BTF 
module 268 to simulate the barb-barbule 
mesostructure. In one embodiment, a quadrilateral 
strip 404 including a plurality of polygons 406 is 
placed along the polyline 4 02 and the texture is 

20 added to each of the polygons 406. The local lighting 
direction 1 and viewing direction v are calculated at 
every vertex x± of polyline 402 using the local 
coordinate frame at Xj.. At each segment edge across 
the barb polyline 402, a ID texture is created by 

25 looking up color values from BTF module 268 generated 
using the feather texture image in window 306 using 
the directions v(xi) and 1 (xi) . Thus, (n+1) ID 
textures are obtained. These textures are combined 
with the texture of the feather in each of the 

30 polygons 406 to render the barb 400 by multi- 
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texturing (application of several textures on the 
same object) and alpha-blending (altering of the 
transparency of an object), as is known in the art. 
When sampling the BTF with a ray tracer, parameters 
can be adjusted so that the BTF gives a "hard" or 
"soft" appearance to the feather. Furthermore, 
occlusions and specularities can be rendered as 
caused by the barb mesostructure . 

Once individual feathers have been 
rendered, the feathers can be placed on a surface. 
FIG. 8 illustrates a framework 500 for placing 
feathers on a surface such as a bird. In framework 
500, several inputs are provided to a bird rendering 
system 502 in order to output a model bird having 
feathers. The inputs include wings and tail skeletons 
504, key feather positions and growing directions 
506, a model of a surface or bird 508 and feather 
geometry 510. Using wings and tail skeletons 504, a 
user will place key feathers and orient growing 
directions of the feathers on the wings and tail 
skeletons. As described below, the wings and tail 
skeletons 504 and the key feather positions and 
growing directions on the wings and tail skeletons 
are provided to an interpolation module 512 within 
bird rendering system 502. 

In order to place contour feathers, key 
feather positions and growing directions are provided 
on the model of a surface or bird 508. Inputs 506, 
508, as well as the feather geometry 510 are provided 
to a recursion module 514. As discussed below, 
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recursion module 514 automatically places feathers 
onto a surface given the feather geometry 510 in 
order to detect and correct collisions between 

adjacent feathers. 

Feathers are placed on the wings and the 
tail using skeletons 504 and key feather positions 
and growing directions 506. The feathers on a wing 
include the primaries, secondaries, humerals, primary 
coverts, and secondary coverts (these feathers are 
rooted on the scapula, ulna/radius, metacarpus and 
phalanx respectively) . As shown in FIG. 9, a polyline 
530 of line segments is used to represent the wing 
skeleton. Similarly a quadrilateral 532 is used as 
the skeleton for the tail. A bird has about 9 to 11 
primaries, 6 to 24 secondaries, and 8 to 24 tail 
feathers. The user can specify the numbers of 
feathers of each type and edit key feathers on the 
wing and key feathers on the tail using polyline 530 
and quadrilateral 532. Interpolation module 512, 
using polyline 530 and quadrilateral 532, generates 
other feathers by interpolation. 

FIG. 9 further illustrates placement of 
feathers on a wing. A model 534 includes a polyline 
536, a surface model 538 and a plurality of segments 
540 indicating where on polyline 536 feathers are to 
be placed. Model 542 illustrates the placement of the 
feathers on the polyline 536 and model 544 
illustrates small wing features placed on surface 
model 538. 



Contour feathers are the feathers that 
cover the body of a bird. Given a polygonal model 508 
describing a bird's body (without feathers), feathers 
of different sizes and shapes are placed on the 
model. The large number of feathers on a bird makes 
it difficult to manually place and edit individual 
feathers. A user can specify a number of key feathers 
and their growing directions using input 506 and 
system 502 the system automatically generates a full 
coverage of the bird based on the key feathers using 
recursion module 514. This full coverage is performed 
by recursion module 514 in three steps: 

a) re-tile the polygonal model 508 to generate 

feather growing positions, 

b) interpolate the key feather growing 
directions to all feather growing points to 
get an initial growing direction at each 
point , and 

c) recursively determine the final feather 

growing direction at every feather growing 
point, with collisions between feathers 
detected and rectified. 

The output of recursion module 514 is a 
feather placement map indicating feather growing 
positions and directions. The feather shape 
parameters (i.e. outline curves 256, rachis curve 258 
and barb curves 260) can be estimated from that of 
nearby key feathers or provided with model feather 
input 510. These shape parameters are used to 
generate a simplified geometry for each feather. This 
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simplified geometry is used for collision detection 
in step (c) above. 

The vertices of the given polygonal model 
508 are difficult to use as feather growing 
5 positions. Feathers at different parts of a bird have 
different sizes, for example small feathers need to 
grow densely in order to cover the bird's skin. In 
addition, feathers tend to distribute evenly in a 
region of constant feather density. Vertices of a 

10 conventional polygonal model often do not have evenly 
distributed vertices or account for increased feather 
density in positions where the feathers are small. To 
address this problem the polygonal model 508 can be 
retiled using a suitable re-tiling algorithm as is 

15 known in the art. This re-tiling creates a polygonal 
model whose vertices are evenly distributed over a 
region of constant density. 

A simple technique may also be used for 
adjusting vertex density based on curvature. This 

2 0 technique can be used to control vertex density based 
on the sizes of feathers. In one embodiment, the 
sizes from the key feathers are estimated using 
Gaussian radial basis functions, where a radius is 
defined as a distance over a surface, as computed 

2 5 using a shortest path algorithm. The user has control 
over the spatial extent and weight of each basis 
function. This estimation scheme is used for creating 
vector fields on a polygonal surface. Ultimately, 
vector fields for vertices are created for the 

30 polygonal model. After re-tiling, the vertices of the 



new polygonal model are the feather growing 
positions . 

From the growing directions of key 
feathers, initial growing directions are calculated 
at all vertices using Gaussian radial basis functions 
as described above. The initial growing directions 
tend to cause inter-penetration of feathers because 
these directions are derived without any 
consideration to the shapes of the feathers. To 
determine the final growing directions, collision 
detection is performed on the feathers based on a 
simplified geometry and the feather growing 
directions are adjusted accordingly. Because of the 
large number of feathers and the complex shape of a 
bird's body, a collision detection between every pair 
of feathers is likely to be very expensive. To 
address this problem two strategies are adopted. 
First, feathers are grown in an orderly fashion 
according to the initial growing directions. Second, 
local collisions between neighboring feathers are 
only considered because collisions rarely happen 
between feathers far away from each other (two 
feathers are neighboring feathers if their growing 
positions are connected by an edge of a polygon) . A 
recursive collision detection algorithm is used in 
recursion module 514 to implement these strategies. 

As FIGS. 10A-10D illustrate, vertices 
around each vertex v are classified into two groups 
according to the initial growing direction o v for v. 
The first group includes vertices in the growing 
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direction and the second group includes vertices not 
in the growing direction. 

After the vertices around every vertex are 
classified, recursion module 514 invokes the 
5 following recursive collision detection algorithm at 

each vertex. 

FindGrowingDirection (v) 

{ 

If the growing direction at v has already 
20 been found, return ; 

For each vertex V in the growing direction 

v, FindGrowingDirection (V ) ; 
While feather (v) collides with feather (V) 

for some vertex in the growing 
15 direction, adjust the growing 

direction at v; 
For each vertex V not in the growing 

direct ion , FindGrowingDirec t ion ( v ' ) ; 

} 

20 Here feather (v) and feather (V) are the 

feathers at vertices v and V, respectively. A vertex 
v on the re-tiled model 508 with initial growing 
direction o v is shown in FIG. 10A. The vertices around 
v are vertices v 0 through v 5 . Among these, v x and v 2 

25 are in the growing direction, whereas v 3 , v 4 and v 5 
are not in the growing direction. According to the 
algorithm above, the final growing directions at v a 
and v 2 is first determined by recursion as illustrated 
at FIG. 10B. Based on the final growing directions at 

30 vi and v 2 as well as the shapes of feathers at v, v x 
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and v 2 , collision between the feathers can be detected 
and the feather growing direction at v can be 
adjusted by rotating feather v toward the surface 
normal at v as shown at FIG. IOC. Rotation can be 
5 performed in small increments, so that rotation can 
stop as soon as there are no collisions detected. The 
growing direction at v is then final and the growing 
directions for v 0/ v 3/ v 4 and v 5 can be processed 
through using the recursion algorithm above as shown 

10 at FIG. 10D. For faster collision detection, a 
simplified geometry can be used for each feather. 

FIG. 11 illustrates a sample bird 550 
having feathers placed according to initial growing 
directions and a sample bird 552 having feathers 

15 rendered according to the algorithm above. In sample 
bird 550, there are numerous inter-penetrating 
feathers and many feathers that grow inside the 
bird's skin as illustrated at 551. When using the 
recursion algorithm above with recursion module 514, 

20 sample bird 552 having a realistic rendering is 
achieved. 

Although the present invention has been 
described with reference to particular embodiments, 
workers skilled in the art will recognize that 
25 changes may be made in form and detail without 
departing from the spirit and scope of the invention. 



